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Understanding  and  control  of  the  transfer  of  oxygen  and  other 
heteroatoms  to  and  from  organic  and  inorganic  species  is  pursued  in 
academic  as  well  as  in  industrial  research.^  Among  commercially  feasible 
heterogeneously  catalyzed  reactions,  the  class  of  selective  oxidations  of 
hydrocarbons  is  very  important.  The  process  of  allylic  ammoxidation  of 
propylene  to  acrylonitrile  on  various  mixed  metal  oxides^  is  an  important 
reaction  and  4.3  million  tons^  of  the  nitrile  are  so  produced  per  year 
worldwide.  The  epoxidation  of  ethylene  to  ethylene  oxide  on  a  silver 
surface^  is  also  very  used  in  the  industry  (5.8  million  tons  per  year^). 
Generally,  for  these  very  important  and  useful  reactions  the  underlying 
mechanisms  of  the  selective  oxidations  are  not  well  understood.  Some 
theoretical  investigations  of  acrolein  production  from  propylene  using  a 
mixed  molybdenum-bismuth-oxide  catalyst^  and  for  the  adsorbtion  of 
ammonia  on  VjOs,  an  important  step  in  the  ammoxidation  reaction,2‘^have 
been  published.  A  study  of  ethylene  epoxidation  on  a  Ag(llO)  surface^  has 
also  apperead  recently.  Perhaps  the  most  intriguing  process  among  the 
heterogenously  catalyzed  selective  oxidation  regions  is  the  synthesis  of 
maleic  anhydride  from  «-butane  over  a  vanadium  phosphorous  oxide 
surface  in  the  presence  of  molecular  oxygen  [equation  (1  )]7 


♦  3.5  O2  (g) 


(V01;P2  07. 


This  reaction,  being  a  14-electron  oxidation  involving  an  abstraction  of 
eight  hydrogen  atoms  and  an  insertion  of  three  oxygen  atoms,  has  been  the 
subject  of  a  large  number  of  reports^'*  and  will  be  the  matter  discussed  in 
this  paper. 

The  solid  state  chemistry  of  the  vanadiu-o  phosphorous  oxide  system  is 
characterized  by  a  large  number  of  ^talline  phases^  that  easily 
interconvert  upon  reduction  and  oxidation.^*' The  active  vanadium 
phosphorous  oxide  phase  in  the  oxidation  of  «-butane  has  been  disputed 
in  the  literature.^^'^'^®  Generally,  it  is  agreed  that  the  presence  of  the 
vanadyl  pyrophosphate  phase,  (V0)2P207,  in  the  catalyst  is  necessary  for  the 
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reaction  to  proceed.^  ^  An  analysis  of  a  commercially  used  vanadium 
phosphorous  oxide  catalyst  showed  only  the  precense  of  one  crystalline 
phase  by  X-ray  diffraction  measurements,  namely  the  vanadyl 
pyrophosphate  phased*  Furthermore,  it  revealed  an  oxidation  state  for 
vanadium  of  4.00-4.03  and  a  P;V  ratio  of  1.00±0.02,^®  consistent  with  the 
assignment  of  (VO)2P207as  the  active  phase.  Most  of  the  mechanistic 
studies  are  thus  performed  using  a  dean  vanadyl  pyrophosphate  surface  as 
the  catalyst.  Here  also  we  will  concentrate  on  this  particular  vanadium 
phosphorous  oxide  structure. 

Vanadyl  pyrophosphate  has  a  layered  structure.  The  layers  consist  of 
pairs  of  edge-sharing  vanadium-oxygen  octahedra  linked  together  in  the 
plane  by  pyrophosphate  units,  1.  The  planes  are  linked  together  through 
the  axial  oxygens  of  the  edge-sharing  octahedra.  A  tetragonal  distortion  of 
the  octahedral  environment  of  vanadium  in  the  double  chain  direction  is 
observed  by  the  X-ray  crystallographic  investigations.^^  The  two  vanadyl 
groups  hereby  formed  are  placed  trans  to  one  another  within  one  pair  of 
edge-sharing  vanadium  octahedra,  2.  The  interlayer  bonding  can  be 
described  as  a  weak  Lewis  acid  ->  Lewis  base  interaction  (V=0  — >V=0).®^ 
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It  is  very  important  to  obtain  knowledge  of  the  surface  chemical 
structure  of  (VO)2P207  in  order  to  trace  a  likely  mechanism  for  the  selective 
oxidation  of  «-butane  to  maleic  anhydride.  Satsuma  et  al.  have  described 
the  surface  structure  in  terms  of  the  number  of  surface  V=0  groups  and 
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surface  acid  sites,  using  various  physical  ana  chemical  methods.^^  They 
were  able  to  show  that  the  vanadyl  group,  V=0,  is  presented  to  a  large 
extent  on  the  surface  of  vanadyl  pyrophosphate.  The  surface  add  sites  were 
identified  as  Brondsted  add  sites,  such  as  P-OH,  whereas  a  small  number  of 
Lewis  acid  sites  were  found  assodated  to  the  surface  vanadivun  atoms.  In 
this  study  we  will  apply  one  layer  of  the  vanadyl  pyrophosphate  struture  as 
a  model  for  the  active  surface.  One  vanadyl  dimer  contains  both  the 
surface  V=0  unit  as  well  as  a  Lewis  add  site,  namely  at  the  vanadyl  group 
pointing  towards  the  bulk. 

The  rate-determining  step  in  the  oxidation  of  butane  on  a  vanadyl 
pyrophosphate  surface  has  been  shown  to  be  the  activation  of  the 
methylene  carbon-hydrogen  bond(s).®®''^  In  order  to  trace  the  steps  of  the 
reaction  that  are  responsible  for  the  large  selectivity  to  maleic  anhydride, 
the  oxidation  of  the  thereby  formed  unsaturated  spedes,  such  as  1 -butene 
and  1,3-butadiene  has  been  very  intensively  studied.^^*^'^^  For  these  species 
the  rate-controlling  step  towards  maleic  anhydride  is  the  insertion  of 
oxygen  into  the  hydrocarbon. This  is  very  different  from  other  olefinic 
selective  oxidations,  as  for  example  the  propylene  to  acrolein  synthesis.^ 

Concerning  the  oxygen  source  for  the  maleic  anhydride  reaction,  it  also 
differs  from  other  selective  oxidation  reactions.  Usually,  the  bulk  of  the 
catalyst  is  involved  in  the  oxidation  of  the  organic  fragment  by  undergoing 
a  reduction  in  the  promotion  of  oxygen  to  the  active  surface.^  ^  In  the 
formation  of  acrolein  from  propylene  on  a  mixed  molybdenum-bismuth- 
oxide  surface  the  availability  of  gas-phase  oxygen  is  not  an  absolute  need.^'^^ 
For  this  reaction  it  has  been  shown  by  labelling  studies  that  bulk  oxygen 
atoms  are  transferred  to  the  hydrocarbon.^^  In  the  case  of  vanadyl 
pyrophosphate  bulk  oxygens  are  not  involved  in  the  reaction^®  and  only 
small  amounts  of  maleic  anhydride  are  formed  in  the  absence  of  02(g).^®'®®'® 
The  main  product  of  1-butene  oxidation  was  thus  crotonaldehyde,  obtained 
by  an  oxy-dehydrogenation  and  an  allylic  oxidation.^®'’ 

How  the  molecular  oxygen  becomes  activated  at  the  vanadyl 
pyrophosphate  surface,  and  in  which  form  (gas-phase,  physisorbed  or 
surface-lattice)  it  becomes  incorporated  into  the  products,  is  still  an 
unanswered  and  disputed  question.^'®  It  has  been  rationalized  that  the 
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vanadyl  pyrophosphate  surface  can  adsorb  one  oxygen  molecule  per  two 
surface  vanadium  atoms.®®  Also,  by  using  butadiene  as  starting  material 
Centi  and  Trifird  were  able  to  show  that  the  furan-oxygen  is  introduced 
through  an  attack  of  an  electrophilic  oxygen  on  the  1,4-positions  of 
butadiene  to  give  2,5-dihydrofuran.'®*  Recently,  Busca  and  Centi®^  detected 
at  least  three  different  intermediates  besides  the  2,5-dihydrofuran,  namely 
an  unsaturated  lactone  and  two  structures  of  adsorbed  furan.  Two 
competing  reaction  paths  could  be  distinguished,  as  shown  in  Scheme  1. 
There  is  postulated  a  fast  path,  involving  an  asymmetric  unsaturated 
lactone  (lower  reaction  path)  and  a  slower  one,  passing  through  a  furan 
intermediate  (upper  reaction  path).  The  slow  path  was  shown  to  be  the  less 
selective  way,  giving  rise  to  the  formation  of  CO  and  CO2.®*’ 


CO  •  CO 


Scheme  1 

In  this  study  we  will  focus  on  the  lower  more  selective  path.  The 
"upstanding"  Ocoordinated  form  of  adsorbed  furan  is  more  stable  by  ~2-3 
eV  than  the  "lying  down"  ^-coordinated  form,  depicted  in  the  upper  path. 
An  analysis  of  the  bonding  situation  between  a  vanadyl  pyrophosphate 
surface  and  the  various  adsorbates  (butadiene,  23-dihydrofuran  and  the 
lactone)  is  made  in  order  to  trace  the  most  likely  adsorption  geometry  of 
these  intermediates.  The  nature  of  the  active  oxygen  is  also  discussed  and  a 
mechanism  for  the  selective  oxygen  transfer  is  proposed.  Throughout  the 
study  the  extended  Hiickel^^  approximation  is  used.  Due  to  limitations  in 
computational  time  a  cluster  model  is  applied  for  the  more  complex 
structures;  otherwise  tight  binding^®  calculations  are  performed.  Further 
details  of  the  calculations  are  given  in  the  Appendix. 
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Let  us  start  with  an  analysis  of  the  clean  vanadyl  pyrophosphate  surface. 
Although  the  structure  of  the  siirface  is  disputed  in  the  litterature^*'^''^  we 
will  here  limit  ourselves  to  the  (020)  plane  of  (VO)2P207  as  the  active 
surface.^^'i^  It  has  been  proposed  that  the  active  sites  are  located  at 
interfaces  between  different  crystal  phases.^  Though  a  crystal  of  one  phase 
only,  the  vanadyl  pyrophosphate  phase  serves  well  as  a  catalyst  for  the 
selective  oxidation  of  m -butane  to  maleic  anhydride.®  Surface  sensitive 
experiments  conclude  that  vanadyl,  V=0,  sites  are  found  at  the  surface  of  a 
(V0)2P207  crystal  and  also  that  surface  Lewis  acid  sites  are  foimd.^3  One  pair 
of  edge-sharing  VO5  units  contains  both  of  these  functional  groups,  3. 


Lewis  acid 
site 


vanadyl  surface  group 


(020)  plane 


0 


b-direction 
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The  surface  of  V2O5  has  been  analyzed  theoretically.^*^  The  structure  of 
V2O5  consists  also  of  tetragonally  distorted  edge-sharing  octahedra.  Because 
of  the  highly  ionic  structures  for  V2O5  and  (VO)2P207  the  bonding  within 
both  structures  is  relatively  localized.  Though  differences  in  the  properties 
of  the  surface  layers  are  computed.^*^'^^ 

Due  to  limitations  in  computational  time  we  will  only  use  one  layer  of 
the  vanadyl  pyrophosphate  structure  as  a  model  for  the  active  surface. 
Usually,  in  the  case  of  metal  single  crystal  surfaces  three  layers  or  more  are 
recommended,2i  but  for  more  ionically  bonded,  therefore  localized,  systems 
a  one  layer  model  has  been  applied  succesfully  as  e.g.  in  the  analysis  of  the 
adsorption  of  organo-rhodium  compounds  on  a  rutile  surface. 22  The  unit 
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cell  within  one  (020)  plane  of  vanadyl  pyrophosphate  is  rather  large,  due  to 
the  alternating  orientations  of  the  pyrophosphate,  groups.  To  avoid 

this  enormous  unit  cell  we  make  a  one^imensional  approximation.  In  the 
following  we  will  compare  characteristic  features  of  two  one-dimensional 
chains  made  from  pairs  of  edge-sharing  VO5  units.  The  dimension  taken 
into  account  is  the  one  along  the  b-direction  of  the  imit  cell,  meaning  along 
the  direction  of  the  edge-sharing  vanadyl  groups.  The  first  model  consists 
of  the  bare  vanadyl  units,  4.  The  second  structure  applied  attaches  PO3 
units  as  models  for  the  pyrophosphates,  5.  The  charge  of  each  structural 
unit  is  chosen  so  as  to  give  vanadium  an  oxidation  state  of  +4  (all  oxygens 
are  in  the  -2  oxidation  state). 
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Small  defects  in  the  (VO)2P207  structure  give  rise  to  a  small  quantity  of 
centers.^*  Some  authors  have  attributed  to  these  the  large  selectivity 
for  maleic  anhydride.7»i5b,23  ^  calculation  with  vanadium  in  the  oxidation 
state  +5  is  thus  performed  in  order  to  investigate  how  the  oxidation  state 
can  influence  the  adsorption  processes  regarded. 

The  distance  between  two  nearest  oxygens  from  two  neighbouring  unit 
cells  is  rather  large  (the  average  oxygen-oxygen  separation  within  a 
pyrophosphate  unit  is  2.50A^^)  so  we  expect  very  little  inter  unit  cell 
interactions.  For  this  reason  we  did  calculations  on  the  parent  compounds 
to  the  two  chains.  A  calculation  with  P(OH)3  units  as  models  for  the 
pyrophosphates,  6,  was  also  carried  out  for  the  +4  oxidation  state  of 
vanadium. 
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Important  results  are  listed  in  Table  1  for  the  calculations  with 
vanadium  in  a  +4  oxidation  state  and  in  Table  2  for  an  oxidation  state  of  +5. 


Table  1 
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For  both  oxidation  states  considered  the  differrent  models  used  gives 
consistent  results.  For  the  set  of  calculations  with  vanadium  in  the  +4 
oxidation  state  all  models  give  a  net  charge  on  vanadium  of  approximately 
+2.0  and  of  -1.26  for  the  double  bound  oxygens.  Also,  the  vanadium- 
vanadium  (V-V)  overlap  population  is  computed  to  be  small  and  positive 
in  all  models.  The  overlap  population  of  the  vanadium-oxygen  double 
bond  is  also  almost  unchanged  from  one  model  to  the  other,  varying  from 
0.758  to  0.771. 

The  other  set  of  calculations  with  two  electrons  less,  giving  vanadium 
in  the  +5  oxidation  state,  also  comes  out  with  consistent  numbers  for  the 
charateristics  listed.  The  charge  on  vanadium  is  about  +3.0  in  all  three 
models.  The  overlap  population  of  the  vanadium-oxygen  double  bond 
remains  unchanged  compared  to  the  calculations  for  the  +4  oxidation  state. 
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Electronic  properties  for  the  different  models  used  for  the  vandyl  pyrophosphate  surfaces  with  vanadium  in  a  +4  oxidation  state. 


Figure  1:  Crystal  Orbital  Overlap  Population  (COOP)  for  ths  V=0  bond  in  4 
(a)  and  5  (b). 
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The  V-V  overlap  population  has  now  settled  on  a  small  negative  value. 

The  consistency  between  the  models  used  is  also  found  in  the  V=0 
COOP2'*  (Crystal  Orbital  Overlap  Population)  curves  for  4  and  5,  as  seen  in 
Figure  1  (a)  and  (b),  respectively.  Both  have  a  bonding  peak  at  ~  -15  eV  and 
two  peaks  contributed  by  orbitals  with  antibonding  V=0  character  above 
the  Fermi  level  at  ~  -9.5  eV  and  -8.0  eV.  The  fine  structure  of  the  peaks  is 
different  because  of  the  changed  environment,  but  the  main  features 
remain  the  same.  We  will  discuss  below  the  origin  of  these  features. 


Figure  1 


The  removal  of  two  electrons  from  the  system  by  going  from  the  +4 
oxidation  state  of  vanadium  to  the  +5  oxidation  state  is  accompanied  by  a 
decrease  in  the  small  V-V  overlap  population  and  an  increase  in  positive 
charge  located  on  each  vanadium.  These  observations  can  be  rationalized 
by  an  analysis  of  some  DOS^^  (Density  of  States)  and  COOP  curves.  Some 
plots  for  structure  5  is  shown  in  Figure  2  to  explain  the  sign  and  magnitude 
of  the  V-V  overlap  population  calculated. 


Figure  2 


The  Fermi  level  of  5  is  located  at  -10.2  eV.  The  plot  shown  in  Figure  2 
(a)  is  the  contribution  from  the  vanadium  d^2  orbitals  to  the  total  DOS. 
We  see  that  -10%  of  the  states  are  found  in  the  -14.5  eV  to  -17.0  eV  range. 
These  states  cause  a  bonding  and  an  antibonding  peak  in  the  V-V  COOP 
curve,  as  depicted  in  Figure  2  (b).  The  antibonding  peak  is  larger  in 
magnitude  than  the  bonding  one,  so  the  total  V-V  overlap  population 
contributed  by  these  lower  peaks  is  negative,  as  calculated  for  the  +5 
oxidation  state  structures.  The  major  part  (90%)  of  the  d^2.y  2  states  is  located 
around  the  Fermi  level.  These  states  also  create  a  bonding  and  an 
anti  bonding  V-V  p)eak  in  the  CCX)P.  The  Fermi  level  cuts  at  an  energy  half 
way  through  this  bonding  peak  explaning  the  small  positive  overlap 
population.  This  indicates  only  a  small  coupling  of  the  d  electron  located  at 
each  vanadium  center  in  agreement  with  the  identification  of  a  triplet 
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ground  state  by  EPR  measurements.®*  By  going  to  the  +5  oxidation  state  of 
vanadium  two  electrons  are  removed  form  each  structural  unit.  A  loss  of 
0.96  electrons  for  each  d^2.y2is  computed.  The  Fermi  level,  located  at  -10.7 
eV,  is  now  far  below  the  upper  V-V  bonding  peak,  justifying  the  decrease  in 
V-V  overlap  population. 

The  appearance  of  two  peaks  in  the  EXDS  for  dj^2.y2  is  a  consequence  of 
mixing  with  the  p  orbitals  of  the  bridging  oxygens.  The  DOS  for  p^^  and  py 
are  plotted  in  Figure  2  (c)  and  (d).  We  see  a  large  contribution  in  the  -14.0 
eV  to  -17.5  eV  range  and  a  smaller  one  (-5%)  around  the  Fermi  level.  A 
schema tical  drawing  of  the  interaction  of  d^2.y2  with  the  p^  and  Py  of  the 
bridging  oxygens  is  shown  in  7.  The  d^2.y  2  orbitals  of  two  vanadiums 
interact  weakly  giving  a  bonding  and  an  antibonding  combination.  Each 
of  these  interact  with  the  p  orbital  of  the  bridging  oxygens  with  the  right 
symmetry.  The  four  combinations  obtained  explain  the  variations 
observed  in  the  V-V  COOP  curve. 
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The  bonding  and  antibonding  combinations  of  the  vanadium  d^y 
orbitals  interact  in  a  similar  way  with  p^  and  Py  of  the  bridging  oxygens. 
The  EXDS  for  d^y  is  plotted  in  Figure  2  (e).  Two  major  peaks  are  found 
above  the  Fermi  level  at  -5.0  eV  and  -6.8  eV.  These  are  seen  in  the  V-V 
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COOP  (Figure  2b)  as  a  bonding,  8,  and  an  antibonding  peak,  9. 

Figure  3  displays  some  important  plots  for  the  understanding  of  the 
vcmadium-oxygen  double  bond.  Figure  1  (b)  showed  the  COOP  of  the  bond 
in  question.  Bonding  states  are  fotmd  from  -14  eV  to  -16  eV  and  the 
antibonding  states  are  located  above  the  Fermi  level  with  two  major  peaks 
at  -9.5  eV  and  -8.0  eV.  The  vanadium  orbitals  involved  in  this  bonding  are 


Figure  3 


depicted  in  Figure  3  (a)  (d^i)  and  in  Figure  3  (b)  (d^^,  dy^).  The  bonding 
states  are  made  up  of  interactions  of  d^^  ^^d  dy^  on  the  metals  with  p^  and 
Py  of  the  oxygens  and  of  with  oxygen  p^  and  s.  The  antibonding  peaks 
are  easily  distinguished  with  help  from  the  EXDS's.  In  Figure  3  (b)  the  -10.0 
eV  peak  just  above  the  Fermi  level  is  the  antibonding  n  type  interaction,  10, 
whereas  the  stronger  o  interaction  has  its  antibonding  peak  at  a  higher 
energy,  -8.0  eV,  11  (Figure  3(a)).  No  vanadium-oxygen  states,  bonding  or 
antibonding,  are  located  right  below  the  Fermi  level  explaining  the 
unaffected  V=0  overlap  population  when  the  oxidation  state  of  vanadium 
is  changed. 

The  results  of  these  one-dimensional  calculations  are  fully  consistent 
with  the  picture  obtained  from  calculations  on  the  parent  clusters.  For  the 
more  complex  structures  where  it  is  harder  to  resolve  the  DOS's  and 
COOP'S  we  will  therefore  use  dusters  as  models  for  the  surface. 
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In  the  following  we  will  focus  on  the  oxygen  transfer  from  the  vanadyl 
pyrophosphate  surface  to  butadiene.  The  hydrogen  activation  and 
abstraction  steps  going  from  n-butane  to  butadiene  will  not  be  considered 
in  this  study.  Experimentally  it  is  believed  that  butadiene  and  n-butane 
pass  through  the  same  intermediates  on  the  way  to  maleic  anhydride,  the 
first  intermediate  being  2,5-dihydrofuran.7a«8 


The  formation  and  bonding  of  2.5-dihvdrofuran. 


It  has  been  proposed  that  the  furan  oxygen  is  introduced  through  an 
attack  of  an  electrophilic  oxygen  at  the  1,4-positions  of  butadiene.  This 
statement  is  in  agreement  with  a  simple  orbital  interaction  diagram 
between  the  V=0  group  of  the  (V0)20g®'  duster  and  butadiene,  as  shown  in 
Figure  4.  The  oxygen-carbon  distance  is  set  to  1.60  A  and  the  1,4-positions  of 
butadiene  are  at  the  same  x-coordinate  as  the  V=0  unit  of  the  surface.  The 
left-hand  side  shows  the  frontier  orbitals  of  (VO)20^®‘.  We  recognize  the 
bonding  and  antibonding  combinations  of  the  vanadium  d^2.y2  orbitals  as 
the  HOMO  and  LUMO,  nearly  degenerate  in  this  molecular  calculation. 
Right  above  the  LUMO  are  located  antibonding  combinations  of  d^  and  dy^ 
with  px  and  Py,  as  in  the  one-dimensional  calculations.  These  V-O 
antibonding  orbitals  are  electrophilic  in  character  when  interacting  with 
the  HOMO  of  butadiene. 


Figure  4 


A  decrease  in  overlap  population  for  the  vanadium-oxygen  double  bond 
from  0.761  to  0.495  is  computed,  verifying  that  a  substantial  fraction  of  V=0 
antibonding  orbitals  becomes  populated  upon  adsorption.  A  decrease  in 
overlap  population  for  the  carbon  carbon  double  bonds  after  adsorption  is 
seen  too,  from  1.465  to  1.137,  due  to  a  depletion  of  0.67  electrons  from  the 
HOMO  and  a  donation  of  0.19  electrons  to  the  LUMO  of  butadiene.  A 
vertically  oriented  butadiene  is  also  suitable  for  interaction  with  the 
vanadyl  oxygen,  but  the  energy  of  the  product  is  disfavored  by  ~3  eV 
compared  to  the  geometry  shown  in  the  center  of  Figure  4.  The  product  of 
this  [4+2]  type  reaction  rearranges  to  the  2,5-dihydrofuran  structure.  We 
have  examined  some  different  orientations  of  2,5-dihydrofuran  (DHF)  in 
order  to  trace  a  likely  orientaton  of  the  intermediate  in  the  maleic 
anhydride  synthesis.  Four  orientations  of  DHF  adsorbed  on  a  (VOljO^®' 
cluster  were  tried.  Two  calculations  with  the  DHF  ring  oriented  parrallel  to 
the  surface,  12,  and  two  calculations  with  the  DHF  ring  perpendicular  to  the 
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surface  plane,  13,  were  carried  out.  From  an  energetic  point  of  view  the  two 
structures  12  are  disfavored  by  7  and  9  eV,  respectively,  compared  to  13,  and 
they  will  not  be  considered  further  in  this  study.  The  two  "upstanding" 
forms  give  similar  energies.  Calculations  for  these  struchires  were  carried 
out  on  the  simple  one-dimensional  surface,  5,  to  reveal  whether  or  not  a 
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preference  in  the  orientation  of  DHF  could  be  found.  Table  3  lists  some 
results  for  the  two  structures  examined,  14  and  15.  Very  similar  numbers 
are  computed,  though  the  structure  with  the  ring  along  the  V-V  direction, 
14,  is  slightly  more  favourable  by  0.3  eV. 


Table  3 


From  these  calculations  it  is  concluded  that  DHF  on  a  vanadyl 
pyrophosphate  surface  takes  an  upright  position,  but  no  specific 
conformation  for  this  orientation  can  be  picked  as  the  most  favourable. 
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Let  us  now  proceed  to  a  description  of  the  bonding  situation  between 
DHF  and  the  vanadyl  pyrophosphate  surface.  For  this  purpose  a  one¬ 
dimensional  chain  of  (VO)206®*  is  applied.  The  distance  between  the 
oxygen  of  DHF  and  vanadium  is  arbitrarily  set  to  2.0  A.^  Figure  5  displays 
some  plots  to  examine  the  V-O  and  the  C-O  bonds  of  16.  Figure  5  (a)  is  the 


Figures 


16 


CCX)P  curve  for  the  V-O  bond.  Compared  to  Figure  1  (a)  especially  the 
bonding  part  has  changed  by  splitting  into  many  separated  peaks.  The 
orbitals  that  make  up  the  five  bonding  V=0  peaks  also  have  either  O-C 
bonding  or  antibonding  character  as  seen  from  the  O-C  COOP  in  Figure 
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Table  3.  Total  energies  and  some  overlap  populations  for  14  and  15. 


5  (b).  The  DOS  curves  (not  shown  here)  reveal  that  the  four  smallest  (in 
amplitude)  occupied  peaks  in  the  two  COOP  curves  are  made  up  from  an 
interaction  of  the  low  lying  V-O  rty  orbitals  with  the  Py  on  the  2,5-carbons  of 
DHF.  The  two  lower  peaks  are  the  bonding  combinations  and  the  upper 
ones  are  the  antibonding  ones  resulting  in  no  net  n  type  bonding  between 
oxygen  and  carbon  in  the  DHF  ring,  17  a.  The  major  peak  around  ~  -15  eV 
consists  of  the  a  type  bonding  orbitals,  17  b.  The  corresponding 
antibonding  orbitals  are  located  above  the  Fermi  level.  Again,  the  orbital 
picture  obtained  from  a  one-dimensional  calculation  is  consistent  with  the 
cluster  calculation.  From  Figure  5  (a)  it  is  seen  that  the  Fermi  level  cuts  at 
an  energy  half  way  through  a  strongly  antibonding  V-O  band.  For  the 
structure  with  two  electrons  less,  vanadium  in  the  +5  oxidation  state,  a 
higher  V-O  overlap  population  is  thus  calculated. 


Bulk  oxygens  are  not  involved  in  the  maleic  anhydride  formation  from 
M-butane7®'®  Only  lattice  oxygens  from  the  very  outermost  layers  are 
transferred. The  activation  and  incorporation  of  molecular  oxygen  into 
the  vanadyl  p5n-ophosphate  surface  is  essential  for  keeping  the  large 
selectivity  for  maleic  anhydride.  Pepera  et  a/.®*  have  showed  that  one 
oxygen  molecule  can  be  adsorbed  per  two  surface  vanadium  centers. 
Within  our  model  this  could  be  one  Oj  molecule  per  pair  of  edge-sharing 
vanadyl  groups.  An  oxygen  molecule  can  coordinate  in  two  ways  to  one 
transition  metal  center,^^  as  either  a  T]*-superoxo  compound,  18,  or  in  the 
Ti^-  coordinated  (peroxo)  way,  19.  The  bending  angle,  0,  of  the  known 
superoxo  complexes  varies  from  -115°  to  -170°.^^  OiUy  structures  with  a 
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peroxo  coordinated  dioxygen  has  been  reported  for  vanadium. Table  4 
contains  some  results  for  and  T[^  coordinated  dioxygen  to  a  chain  of 
(V0)20^®'  units,  20.  Three  bending  angles  of  the  superoxo  structure  were 
used. 


20 
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Table  4 


The  separation  of  vanadium  and  one  oxygen  of  the  dioxygen  fragment 
was  taken  to  be  1.90  A  and  the  oxygen-oxygen  distance  to  1.40  hP  A  very 
flat  minimum  exists  for  the  superoxo  structure  as  a  function  of  the  bending 
angle,  0.  The  total  energy  is  unchanged  from  180°  to  ~150°.  An  angle  of 
90°  is  unfavorable  by  -3  eV.  The  total  energy  of  the  peroxo  form  is  only  0.5 
eV  higher  than  the  superoxo  form.  An  angle  of  150°  for  18  is  applied  in  the 
rest  of  the  study. 

For  both  coordination  modes  a  decrease  in  the  oxygen-oxygen  overlap 
population  is  calculated  upon  adsorption,  from  0.570  in  Ojfg)  to  0.395  and 
0.413  for  the  and  coordinated  species,  respectively.  The  vanadium  to 
oxygen  overlap  population  is  highest  in  the  peroxo  structure. 

The  bond  obtained  between  a  t|2  coordinated  dioxygen  and 
can  be  described  as  an  interaction  of  the  partly  occupied  n*  of  O2  and  the 
unoccupied  antibonding  vanadium  oxygen  d^^  +  p^  orbital,  21.  A  gain  of 
0.12  electrons  is  calculated  for  d^^-  superoxo  structure  the  a  bond  of 

O2  is  involved  in  bonding  too.  A  loss  of  0.12  electrons  from  the  a  orbital  is 
computed.  This  electron  density  is  mainly  transferred  to  622  (0.06  electrons) 
and  P2  (0.08  electrons)  of  vanadium.  Again,  the  orbital  of  Oj  interacts 
with  the  antibonding  d^^  oxygen  p^  orbital,  22. 
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Coadsorption  of  02(g>  and  DHF  on  the  vanadyl  pyrophosphate  surface. 
A  model  for  the  formation  of  maleic  anhydride. 


The  presence  of  two  vanadium  atoms,  located  3.19A  from  each  other, 
on  the  vanadyl  pyrophosphate  surface  leads  to  a  possible  activation  of  02(g) 
at  one  of  the  vanadium  atoms,  while  DHF  remains  adsorbed  at  the  other 
vanadium  atom.  The  are  four  main  combinations  of  a  ti^-superoxo 
adsorbed  molecular  oxygen  and  DHF  as  outlined  in  23  -  26. 


For  the  coadsorption  of  ■n2-p>eroxo  molecular  oxygen  and  DHF  another  four 
main  possible  orientations  of  O2  and  DHF  emerge,  as  shown  in  27  -  30. 

Table  5  gives  the  total  energy  and  some  overlap  populations  of  interest 
for  the  present  study. 


Table  5 
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Table  4.  Total  energies  and  some  overlap  populations  for  18  and  19. 


0 

29 


0 — 0  9 

o.A  /  ^0.  I 

0^  I  ^0^ 


30 


,0  !®' 
0 


0 


It  appears  from  Table  5  that  the  energy  difference  is  only  0.21  eV  between 
the  most  stable  structure  with  Ti‘-superoxo  molecular  oxygen  coadsorbed 
with  DHF  on  the  vanadyl  pyrophosphate  surface  and  the  lowest  energy 
structure  of  ii^-peroxo  molecular  oxygen  and  DHF  on  the  surface.  This 
indicates  that  within  the  present  calculational  method  it  is  not  possible  to 
distinguish  between  these  two  lowest  energy  structures  of  adsorbed 
molecular  oxygen  and  DHF.  In  two  other  pairs  of  geometries  of  molecular 
oxygen  and  DHF  coadsorbed  (24  and  28,  and  25  and  29)  the  difference  in 
energy  is  also  very  small,  whereas  26  is  considerably  higher  in  energy  than 
30.  The  latter  case  is  due  to  electronic/steric  repulsion  between  the 
terminal  oxygen  and  the  two  hydrogens  in  the  2-position  of  DHF.  The 
different  overlap  populations  are  remarkably  constant  for  the  different 
orientations  of  either  a  T|*-superoxo  or  a  Ti^-peroxo  molecular  oxygen 
coadsorbed  with  DHF.  The  only  two  overlap  populations  that  are  different 
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Table  5.  Total  energies,  some  overlap  populations  of  23  -  30  and  the  shortest  O-H  (H  of  DHF)  distance. 


for  the  various  orientations  are  O-H  or  C-H  in  26  and  30;  a  negative  O-H 
overlap  population  is  found  in  26  and  a  reduction  in  the  C-H  overlap 
population  is  observed  for  30. 

We  will  now  proceed  to  a  description  of  how  maleic  anhydride  can  be 
formed  from  coadsorbed  molecular  oxygen  and  DHF  on  the  vanadyl 
pyrophosphate  surface.  We  are  aware  of  the  fact  that  the  quantitative 
results  obtained  by  the  extended  Huckel  calcualtions  cannot  be  sufficiently 
trusted  in  studies  where  several  bonds  are  simultaneously  formed  or 
cleaved.  But,  as  our  intention,  here  in  this  part  of  the  paper,  is  to  give  a 
qualitative  description  of  the  reaction  path,  this  method  should  be  able  to 
provide  us  with  some  useful  trends.  As  a  starting  point  for  our  analysis  we 
will  study  the  rotation  of  DHF  coadsorbed  with  either  a  superoxo  or  a 
peroxo  molecular  oxygen  around  the  axes  outlined  in  a  topview  in  31  and 
32. 


31  32 

When  rotating  one  of  the  two  coadsorbed  molecules  on  the  vanadyl 
pyrophosphate  surface  the  terminal  oxygen  of  a  q^-superoxo  or  one  of  the 
oxygens  in  a  Tj^-peroxo  adsorbed  molecular  oxygen  is  brought  close  to  the 
hydrogen  of  the  C-H  bond  in  the  2-position  at  DHF.  These  approaches  of 
oxygen  to  hydrogen  leads  to  several  electronic  changes. 

In  the  rotation  DHF  starts  out  being  oriented  perpendicular  to  the  V-V 
direction  of  the  surface  whereas  the  oxygen  fragment,  either  superoxo  or 
peroxo,  stays  fixed  along  the  V-V  direction  (31  and  32).  Figure  6  shows  the 
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changes  in  energy  for  the  two  different  systems  as  a  function  of  the  O-H 
distance.  The  rotations  show-n  correspond  to  rotation  angles  from  0°  to  75° 
round  the  V-O  (DHF)  axis. 


Figure  6 


From  Figure  6  it  is  seen  that  in  the  course  of  the  rotation  of  DHF,  when 
coadsorbed  with  a  peroxo  species,  the  approach  of  the  hydrogen  in  DHF 
leads  to  an  increase  in  energy.  It  costs  about  4.6  eV  to  rotate.  In  the 
superoxo  case  about  6.8  eV  is  required  for  the  approach  of  the  hydrogen  to 
oxygen  to  the  same  O-H  distance  (1.15  A),  which  indicates  that  the  former  is 
the  least  unfavorable.  Besides  the  changes  in  energy  several  changes  in 
bonding  is  also  observed:  For  the  rotation  31  the  C-H  overlap  population  in 
DHF  is  reduced  from  0.826  to  0.525,  whereas  the  O-H  overlap  population 
increases  from  0.000  to  0.209  over  the  range  the  O-H  bond  distances  in 
Figure  6.  This  tells  us  that  the  C-H  bond  is  weakened  and  that  a  O-H  bond 
is  made.  The  V-O  (the  oxygen  which  interacts  with  the  hydrogen)  overlap 
population  in  32  is  also  affected  by  the  rotation;  it  decreases  from  0.281  to 
0.234.  In  rotation  32  the  C-H  overlap  population  of  DHF  is  also  reduced 
from  0.826  to  0.286,  whereas  the  O-H  overlap  population  again  shows  an 
increase  from  0.000  to  0.408  for  the  O-H  bond  distance  changes  of  Figure  6. 
This  also  indicates  a  weakening  of  the  C-H  bond  and  an  O-H  bond  being 
formed  too. 

The  change  in  nature  of  the  C-H  bond  of  DHF  and  of  the  interaction 
between  an  oxygen  from  the  peroxo  fragment  and  hydrogen  in  DHF  can  be 
traced  to  a  controlling  orbital  interaction;  the  o^-.h  donates  electron  density 
into  the  cj*o-o  orbital  which  is  the  LUMO  of  the  peroxo  group  that  is 
coordinated  to  vanadium  as  shown  in  33.  The  orbital  donates  0.3 
electrons  into  the  o*q_q  orbital  by  this  interaction. 

The  donation  of  electron  density  from  the  orbital  into  the  o‘o-o 
orbital  thus  leads  to  the  weakening /cleavage  of  the  C-H  bond  and  the 
formation  of  an  O-H  bond.  A  similar  type  of  C-H  activation  has  also  been 
found  in  the  reactions  of  several  types  of  transition-metal  peroxo 
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0-H  (A) 


Figure  6:  Variation  in  the  total  energy  of  31  and  32  when  rotating  around 
the  V-O  (DHF)  axis.  The  O-H  distance  is  measured  from  the  hydrogen  in 
the  2-position  of  DHF  to  the  nearest  oxygen  located  at  the  adjacent 
vanadium  atom. 


complexes;^®  in  the  case  of  transition-metal  peroxo  complexes  the  reaction 
with  e.g.  benzene  leads  to  a  hydroxylation  reaction.^® 
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33 

The  first  step  of  this  proposed  mechanism  causes  a  hydrogen  transfer 
from  the  2-position  in  DHF  to,  probably,  one  of  the  oxygens  of  a  coadsorbed 
Ti  ^-peroxo  oxygen  fragment,  although  we  can  not,  on  the  basis  of  these 
calculations,  exclude  the  Ti^-sup)eroxo  species  being  active  in  the  process. 
When  an  oxygen  from  the  T]^-peroxo  species  interacts  with  the  hydrogen 
the  observed  decrease  in  overlap  population  for  the  V-O  bond  (O  from 
peroxo)  might  indicate  the  this  bond  will  break  by  the  hydrogen  transfer 
reaction,  leading  to  a  t]  *  coordinated  O-O-H  as  shown  in  34.  The  same 
intermediate  is  also  formed  in  the  case  of  an  interaction  between  a  t|’- 
superoxo  dioxygen  coadsorbed  with  DHF. 
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Calculation  of  the  preferred  orientation  of  the  O-H  group  in  34  shows 
that  an  orientation  with  the  hydrogen  pointing  in  the  z-direction  is  0.24  eV 
more  favorable  than  the  one  where  the  hydrogen  points  in  the  opposite 
direction.  With  34  as  intermediate  on  the  reaction  path  the  formation  of 
the  2-hydroxy  derivative  of  DHF  is  a  feasible  process.  The  orientation  of 
DHF,  for  which  the  hydrogen  is  transferred  to  the  oxygen,  is  also  the 
optimum  position  for  interaction  between  O*  and  of  34.  In  this  position 
the  C^-O^  overlap  population  is  0.635  (C^-O^  distance  1.60  A)  showing  that 
significant  C-O  bonding  is  present.  If  the  DHF  derivative  in  34  is  rotated 
away  from  this  optimum  position  (  with  the  DHF  derivative  rotated  15° 
away  from  the  V-V  direction)  to  e.g.  a  position  where  the  DHF  derivative 
is  oriented  perpendicular  to  this  direction,  the  energy  required  for  bringing 
it  back  to  the  optimun  position  is  less  that  2  eV.  Besides  the  formation  of  a 
C-O  bond  the  0^-0^  overlap  population  is  also  reduced  slightly  when  C^  and 
interacts.  With  the  DHF  derivative  oriented  perpendicular  to  the  V-V 
direction  the  0^-0^  overlap  population  is  0.420,  and  0.401  in  the  position 
where  the  optimum  interaction  between  C’-O’  occurs.  The  process 
described  here  gives  then  the  2-hydroxy  derivative  of  DHF  (2FIDHF)  and  an 
adjacent  V-O  function  on  the  surface,  35. 
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The  above  description  of  the  formation  of  35  is  a  two-step  process,  but 
we  can  not  exclude  a  one-step  process,  an  "attachment-insertion"  reaction 
similar  to  the  one  suggested  for  the  hydroxylation  of  C-H  bonds  by  oxene- 
like  reagents.^^ 

The  presence  of  the  V-O  function  adjacent  to  2HDHF  is  also  important 
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for  the  further  reaction.  The  V-O  function  is  located  in  an  optimum 
position  for  interaction  with  the  hydroxy  function  in  2HDHF,  leading  to  the 
asymmetric  unsaturated  lactone.  The  electronic  structure  of  the  V-O 
function  can  be  established  from  the  frontier  orbitals  of  35  shown  in 
Figure  7. 


Figure? 


Figure  7  shows  that  the  V-O  function  contributes  to  the  frontier  orbitals 
with  low  lying  unoccupied  as  well  as  with  relatively  high  lying  occupied 
molecular  orbitals  of  both  p^^  Py  and  p^  character.  These  orbitals  then  have 
the  possiblity  of  interacting  with  either  the  O-H  or  C-H  bond  in  2HDHF  by 
accepting  or  donating  electron  density.  The  interaction  of  the  O-H  and  the 
C-H  bond  in  2HDHF  is  analyzed  as  a  function  of  the  rotation  angle  of 
2HDFIF,  going  from  -30®  to  +30°,  as  shown  in  36,  in  an  attempt  to  elucidate 
the  next  step  in  the  formation  of  the  asymmetric  unsaturated  lactone. 
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In  the  rotation  sketched  in  36,  the  shortest  H^O^  distance  is  1.17  A  and 
the  shortest  H^-O^  is  1.61  A,  indicating  that  the  best  interaction  between 
2HDHF  and  the  oxygen  on  the  adjacent  vanadium  atom  is  the  interaction 
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Figure  7:  The  frontier  orbitals  of  the  vanadium-oxygen  bond  in  35. 


of  with  the  of  the  hydroxyl  group  of  2HDHF.  The  total  energy  of  the 
system  rises  by  about  2  eV  upon  a  rotation  of  +30°.  Several  electronic 
changes  takes  place  upon  rotation:  The  H^-O^  distance  is  1.78  A  in  the 
starting  position  with  the  following  overlap  populations:  C^-O^  0.601,  O^-H^ 
0.603  and  H^-O^  -0.02;  in  the  position  with  the  shortest  distance  (1.17 

A)  these  overlap  populations  have  changed  to:  0.633,  0.366  and  0.154, 
respectively.  These  results  show  that  when  the  C^-O^  bond  increases  in 
strength,  meaning  that  the  C=0  bond  is  beginning  to  form,  the  O^-H^  bond 
is  weakened  and  a  H^-O^  bond  is  being  formed.  The  formation  of  H^-O^ 
bond  can  be  traced  to  interaction  between  the  second  HOMO  shown  in 
Figure  7  and  the  in  the  hydroxyl  group  in  2HDFIF. 

By  the  reaction  path  presented  so  far  we  are  now  near  the  end  of  the 
reaction,  for  we  have  reached  a  hydroxy  group  at  the  adjacent  vanadium 
atom  and  a  2HDHF  radical-like  species  coordinated  to  to  the  other 
vanadium  atom,  as  shown  in  37. 
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With  an  orientation  of  the  system  as  depicted  in  37,  the  last  step  in  the 
formation  of  the  asymmetric  unsaturated  lactone  is  an  abstraction  and  a 
transfer  of  to  O^,  and  a  rehybridization  of  from  sp^  to  sp^.  The  H^-O^ 
distance  in  37,  with  the  bond  lengths  and  angles  used  here  for  the  adsorbed 
species  (see  Appendix)  is  calculated  to  be  1.61  A,  and  for  this  distance  a 
small  H^-O^  overlap  population  is  observed  (0.06),  indicating  that  a 
hydrogen  transfer  might  be  possible.  This  transfer  of  to  O^,  under 
formation  of  H2O  and  the  asymmetric  unsaturated  lactone,  38,  leads  to  a 
decrease  in  the  total  energy  of  the  system  by  0.87  eV,  also  showing  that  the 
process  seems  possible. 
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38 

By  the  process  described  so  far  DHF  has  by  use  of  one  molecular  oxygen 
been  oxidized  to  the  asymmetric  unsaturated  lactone  and  a  water  molecule 
on  the  vanadyl  pyrophosphate  surface.  The  next  step  in  the  formation  of 
maleic  anhydride  on  the  vanadyl  pyrophosphate  surface  is  an  exchange  of 
the  water  molecule  with  molecular  oxygen  in  38  and  a  rotation  of  the 
asymmetric  unsaturated  lactone  by  180°.  This  would  lead  another  CHj 
group  to  be  near  the  newly  adsorbed  and  activated  molecular  dioxygen  at 
the  adjacent  vanadium  atom.  The  formation  of  maleic  anhydride  can  thus 
take  place  in  a  way  similar  to  that  described  above  for  the  formation  the 
asymmetric  unsaturated  lactone.  The  present  results  seem  to  support  the 
observation  by  Busca  and  Centi,  that  the  lower  reaction  path  of  Scheme  1 
represents  a  likely  path  for  the  formation  of  maleic  anhydride  on  a  vanadyl 
pyrophosphate  surface. 
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In  this  study  we  have  analysed  the  electronic  structure  of  the  vanadyl 
pyrophosphate  surface  trying  to  explain  some  of  the  steps  involved  in  the 
transformation  of  n -butane  to  maleic  anhydride  catalyzed  by  this  surface 
using  molecular  oxygen  as  the  oxygen  donor.  First,  an  electrophilic  oxygen 
from  a  surface  vanadyl  group  can  insert  in  the  1,4-positions  of  butadiene  in 
a  concerted  [2+4]  like  cycloaddition  process.  The  thereby  formed  2,5- 
dihydrofuran  fragment  rearranges  itself  to  a  "vertical"  position  where  it  can 
interact  with  a  neighbouring  adsorbed  dioxygen  species. 

The  nature  of  this  dioxygen  fragment  when  adsorbed  on  a  clean  vanadyl 
pyrophosphate  surface  was  analyzed.  The  two  coordination  modes,  (0  = 
150°)  and  r(^,  gave  quite  similar  results,  so  a  definite  assignment  of  the 
geometry  could  not  be  made.  When  interacting  with  a  coadsorbed  DHF 
molecule  the  q^-peroxo  form  seems  to  provide  the  most  likely  structure  of 
the  dioxygen  part,  but  the  q^-superoxo  possibility  cannot  be  excluded.  Both 
dioxygen  structures  upon  abstraction  of  an  2-hydrogen  from  DHF  lead  to 
an  OO-H  surface  species.  This  species  can  transfer  an  OH  group  to  the  DHF 
radical  derivative  forming  a  new  vanadyl  group.  The  newly  formed 
hydroxylated  DHF  can  eliminate  water  upon  reaction  with  the  new 
neighbouring  vanadyl  oxygen  leaving  back,  on  the  surface,  the  asymmetric 
unsaturated  lactone.  Another  molecule  of  O2  can  undergo  the  exactly  same 
reaction  sequence  with  the  5-position  of  the  lactone.  More  detailed 
mechanistic  studies  of  these  later  steps  of  the  proposed  reaction  path  will 
definitely  be  of  great  interest. 
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Table  6:  Atomic  parameters. 


The  extended  Hiickel^^  molecular  and  tight  binding^®  approach  is 
applied  throughout  this  study.  Bond  distances  and  angles  for  the  different 
surface  models  are  taken  as  average  of  the  experimental  values  reported  for 
the  vanadyl  pyrophosphate  structure.^^  This  leads  to  a  V-V  distance  of 
3.188A  and  a  V=0  distance  of  1.628A.  The  geometry  of  2^dihydrofuran  is 
in  accordance  with  the  experimentally  reported  structure.^  Butadiene  has 
C-C  separations  of  1.48 A  and  1.34A  and  all  angles  are  120°.^^  The  dioxygen 
species  have  an  0-0  distance  of  1.40A.^^  The  distance  between  vanadium 
and  an  oxygen  from  the  adsorbed  dioxygen  is  1.90A.^  Atomic  parameters 
are  listed  in  Table  6.  In  the  tight  binding  calculations  properties  averaged 
over  the  Brillouin  zone  were  estimated  by  the  use  of  a  10  fc-point  set 
according  to  the  geometrical  method  by  Ramirez  and  Bohm.^^ 


Table  6 
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